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. . INTRODUCTION 

O I Statistical properties of eigenfunction amplitudes in disordered and chaotic systems 

have attracted considerable research interest recently. Fluctuations of the wave function 

^ ■ amplitudes are believed to determine statistical properties of conductance peaks in 

^ . quantum dots in the Coulomb blockade regime ^' ^' ^' ^' ^, and can be directly measured 

^^ i in the microwave cavity experiments ^' ^. On the theoretical side, the recent progress 

is based on application of the supersymmetry method to the problem of eigenfunction 

CN ' statistics ^' ^. It was found that in so-called zero- mode approximation the distribution 



of eigenfunction amplitudes is correctly described by formulas of the random matrix 



G^ . theory (RMT). Deviations from the RMT predictions were studied in Refs.8, 9, 10, 11. 

'j^ , The present article addresses a problem of correlations of eigenfunction amplitudes. 

A description of correlations in amplitudes of a wave function in relatively close points 
. _ of a chaotic billiard was proposed by Berry ^^ within a RMT-like assumption that the 

pH ■ wave function is a superposition of plane waves with random coefficients. More recently, 

P ■ these correlations were considered in a disordered system within the zero-mode approx- 

imation ^^; the result was later shown ^^ to be equivalent to that of Ref.l2. This result 
is valid for small separations of the two points (less than the mean free path I) . Here 
rS \ we consider such correlations for arbitrary distances. These correlations determine, in 

c^ ■ particular, fluctuations of matrix elements of the (Coulomb) interaction, which are in 

turn important for statistical properties of spectra of quantum dots. Another topic 
addressed in the present paper is that of correlation of amplitudes of different eigen- 
functions. Such correlations, while absent in RMT, appear in a disordered system. 
They become especially strong near the Anderson metal-insulator transition, where 
they play a crucial role in supporting the RMT-like level repulsion. 

This article is based on recent works done in collaboration with Ya. M. Blanter, 
Y. V. Fyodorov, and B. A. Muzykantskii ^' ^^' ^^' ^^' ^^' ^^. 
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EIGENFUNCTION CORRELATIONS IN THE WEAK LOCALIZATION 
REGIME 

In this section, we study the correlations of eigenfunctions in the regime of a 
good conductor ^' ^^' ^^. The correlation function of amplitudes of one and the same 
eigenfunction with energy E can be formally defined as follows: 



a(ri, r2, E) = {\Mri)Mr2)\') ^ ^ A {Y^\Mr,)Mr2)\'S{E - e^) 



(1) 



where ipki^^) and e^ are eigenfunctions and eigenvalues of the Hamiltonian in a partic- 
ular disorder configuration U{r), the angular brackets (. . .) denote averaging over the 
disorder potential f/(r), and A = {J2k^{E — 6^))"^ is the mean level spacing. In the 
case of a real or numerical experiment, calculation of the correlation function (0) would 
include averaging over all states in a relatively narrow energy window around E. An 
analogous correlation function for two different eigenfunctions is defined as 



a{ri,r2,E,uj) 



\Mri)Mr2)\'' 



E,io 



= A^i^-'H ( E \Mri)Mr2)\'S{E - ek)5{E + uj - ei) 

\k^l 

where R2{uj) denotes the two- level correlation function, 

R2iuj) = A^ ( E 6{E - tk)5{E + uo-e{)). 

\ k,l I 



(2) 



(3) 



Eq. (0) defines an overlap of the eigenfunctions ipk and ipi provided they have energies 
close to E with the energy difference equal to u. 

To evaluate a;(ri,r2,-E) and cr(ri, r2, -E, tu) (Ref.l6), we employ an identity 
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A~ia(ri, ra, E)5{uj) + A-^R2{uj)a{ri, ra, E, uj) 

G^ln, n, E)G^{v2. rs, E + u)- G^tn, n, ^)G'^(r2, ra, E 



(4) 



UJ 



where G^'^{y,y' ,E) are retarded and advanced Green's functions and R2{jjj) is non- 
singular part of the level-level correlation function: R2{uj) = R2{^) + S{uj/A). A 
natural question, which arises at this point, is whether the r.h.s. of Eq.(P cannot be 
simply found within the diffuson-Cooperon perturbation theory ^°. Such a calculation 
would, however, be justified only for u ^ A (more precisely, one has to introduce an 
imaginary part of frequency: uj -^ u + iT, and require that F ^ A). Therefore, it 
would only allow to find a smooth in uj part of o"(ri, r2, E, u) for u ^ A. Evaluation 
of a(ri,r2, -E"), as well as of a{ri,r2,E,u) at cu ~ A cannot be done within such a 
calculation. For this reason, we employ a non-perturbative supersymmetry approach 
below. 

The r.h.s. of Eq.(^ can be expressed in terms of the supermatrix cr-model ^^, 
yielding: 



2it' 



A-^a{ri, r2, E)5{uj) + A-V(ri, r2, E, u)R2{u 



= (Trz/)^ [l - Re{Qll{r^)QlUr2))F - Uv, - V2)Re{Q'i{v,)Qll{v^)) ^ 
where kdiv) = (7rz/)^^(ImG''^(r))^ is a short-range function explicitly given by 



kd{r) = exp(-r//) 



(ppr) 



'^ sin^ Ppr, 
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(5) 
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Here (. . .)f denotes the averaging with the action of the supermatrix sigma-model -F[(5]: 

{...)p = lDQ{...)exp{-F[Q]), 

F[Q] = -^Jdr Str[D(Vg)2 + 2t{uj + tO)AQ], (7) 

where D is the diffusion coefficient, u is the density of states, Q = T^^AT is a 
4x4 supermatrix, A = diag(l, 1, — 1, — 1), and T belongs to the supercoset space 
f/(f, f |2)/f/(f |f) X U{1\1). The symbol Str denotes the supertrace defined as StrS = 
Bll — Bj^j- + B^^ — Bjj. The upper matrix indices correspond to the retarded- advanced 
decomposition, while the lower indices denote the boson-fermion one. The action (|^) 
is written for the case of so-called unitary ensemble (broken time reversal symmetry), 
which we consider below. Generalization to a system with time reversal symmetry 
(orthogonal ensemble) is straightforward, and the results are presented in the end of 
the section. Evaluating the a-model correlation functions in the r.h.s. of Eq.(^) and 
separating the result into the singular the singular (proportional to S{lj)) and regular 
at tu = parts, one can obtain the correlation functions a(ri, r2, E) and cr(ri, r2, E, u). 
The two-level correlation function R2{u!) entering Eq.(^ is given in this formalism by 

^2M = ^Re|rfrirfr2 [l - {Qllir,)Qllir,))F 

In the metallic (weak localization) regime, the sigma-model correlation functions 
{Qll{ri)Ql'^{r2))F and {Ql'l{ri)Qll{r2))F can be calculated for relatively low frequencies 
uj <^ Ec with the use of a general method developed in Refs.22, 9, which allows one to 
take into account spatial variations of the field Q. The results are obtained in form of 
expansions in g~^, where g is the dimensionless conductance. First, we restrict ourselves 
to the terms of order g~^. Then, the result for the first correlation function reads as 

(«;(rO«(r.)). = -1 - a ^-P'7'f-- - ^n(r., ..) , (9) 

where s = uj/A + iO. Here the diffusion propagator H is the solution to the diffusion 
equation 

- L'V'H(ri, rs) = ini^)-'[6{T^ - V2) - V-'] (10) 

with the Neumann boundary condition (normal derivative equal to zero at the sample 
boundary), which can be presented in the form 

H(ri, r2) = (nu)-' Y.iDq')-'M^i)M^2), (11) 

q 

with 0q being the eigenfunction of the diffusion operator corresponding to the eigenvalue 
Dq^, q 7^ 0. The first two terms in Eq. (^ represent the result of the zero-mode 
approximation ^^, which takes into account only the spatially constant configurations 
of the field Q{r), so that the functional integral over DQ{r) is reduced to an integral 
over a single matrix Q. The last term is the correction of order g~^. An analogous 
calculation for the second correlator yields: 



{Qlt{ri)Qtl{r2))F 




expuvrs smvrs 
1 +z— ^ ^ 
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n(ri,r2) . (12) 



Now, separating regular and singular parts in r.h.s. of Eq. (^, we obtain the following 
result for the autocorrelations of the same eigenfunction: 

V^{\Mri)Mr2)f)E - 1 = Ur)[l + H(ri, n)] + H(ri, r2), (13) 



and for the correlation of amplitudes of two different eigenf unctions 

V^2(|^,(ri)V^Kr2)r)E,.-l = A;rf(r)n(ri,ri), k^l (14) 

In particular, for ri = r2 we have 

V'{\Mr)tPir)\^)E,u. -1 = 6m + {1 + SmMt, r). (15) 

Note that the result (|13|) for ri = r2 is the inverse participation ratio calculated in Ref. 



9; on the other hand, neglecting the terms with the diffusion propagator (i.e. making 
the zero-mode approximation), we reproduce the result of Refs.l2, 13, 14. 

Eqs. (p!^, (p!5|) show that the correlations between different eigenf unctions are 
relatively small in the weak disorder regime. Indeed, they are proportional to the small 
parameter n(r, r), which is equal in the case of 2D geometry to {L is the size of the 
system) 

Uiv,r) = ing)-HnL/l, 2D, (16) 

with g = 2'KvD. For a quasi-lD wire or strip of the length L, 



2 

n(r, r) = - 
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l^Mi 



< r < L , (17) 



where g = 27ii'D/L, and -82(2;) = x^ — x + 1/6 is the Bernoulli polynomial.^ The corre- 
lations are enhanced by disorder; when the system approaches the strong localization 
regime, the relative magnitude of correlations, n(r, r) ceases to be small. The correla- 
tions near the Anderson localization transition will be discussed in the next section of 
the paper. 

Another correlation function, generally used for the calculation of the linear re- 
sponse of the system, 

7(ri,r2,E,u;) = (^fc(ri)^;(ri)^fc(r2)^;(r2))£;^^ 

= A'R,\u;) (Y^rk{rim^i)M^2)i^:i^2nE - ek)6{E + 00 - e^)) (18) 

can be calculated in a similar way; the result reads 

\/2(^*(ri)^Kri)^fc(r2)^r(r2))i,,. = ^r) + n(ri, ra), k^l. (19) 

As is seen from Eqs. (0), (0)) 0T9|), in the 1/g order the correlation functions 
a{ji, r2, E) and 7(ri, r2, E, uj) survive for the large separation between the points, r ^ I, 
while o"(ri, r2, -E,ci;) decays exponentially for the distances larger than the mean free 
path /. This is, however, an artifact of the g^^ approximation, and the investigation 
of the corresponding tails requires the extension of the above calculation to the terms 
proportional to g"'^. We find that the correlator {Qll{ri)Ql1{r2))F gets the following 
correction: 

HQll{ri)QlUr2))F = -h + 2/4 + exp(2zvr.)/3 - 2i^-^^^^^^^{h - h) 

TVS 

_exp|^ (20) 

2 TTsr 



^ In the 3D geometry, the sum over the momenta q in Eq.(|ll|) determining n(r, r) diverges at large 
q and is determined by the upper cut-ofF, q ~ 1/Z, yielding n(r, r) r^ g~^L/l. This reflects the fact 
that in 3D geometry the truly local (ri = r2) correlations may not be given correctly by the diffusion 
approximation and can depend on microscopic structure of the random potential. For this reason, we 
do not consider local correlations in 3D geometry here. Note, however, that this concerns the global 
geometry of the sample; locally the system can be either of 2D or 3D nature, which determines the 
form of the function kd{r) (e.g, a wire is locally 3D, but has a quasi-lD geometry). 



Here we defined the functions 



/i(ri,r2) = n2(ri,r2), 

/2(ri,r2) = {2V)-' J dT[u\T,ri) + U^{r,r2) 

fs = V-'Jdrdr'U\r,r'), 

U{r^,r2) = V-' J drU{r,r,)U{r,r2). 



(21) 



Consequently, we obtain tlie following results for the correlations of different {k ^ I) 
eigenfunctions at r > /: 



V^\M^^)Mr2)\')E,.. 



l=o(/l-/3-2/4) 



+ 2(/2 - /3 
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(22) 



As it should be expected, the double integral over the both coordinates of this corre- 
lation function is equal to zero. This property is just the normalization condition and 
should hold in arbitrary order of expansion in g~^. 

The quantities /2, /s, and /4 are proportional to g'"^, with some (geometry- 
dependent) prefactors of order unity. On the other hand, /i in 2D and 3D geometry 
depends essentially on the distance r = |ri — r2|. In particular, for / -C r <^ L we find 



/i(ri,r2)=n2(ri,r2) 



(tt^)' 



In^- 



(47r2z/Dr)2 



2D 
3D 



Thus, for I < r '^ L, the contributions proportional to /i dominate in Eg . (|22|) , and we 
get 

V^{\Mri)Mr2)\')E,u.-l = \u\r,,r2) , k^l. (23) 

On the other hand, for the case of quasi- ID geometry (as well as in 2D and 3D for 
r ~ L), all quantities /i, /2, /s, and /4 are of order of l/g"^. Thus, the correlator 
cr(ri,r2, -E,^^) acquires a non-trivial (oscillatory) frequency dependence on a scale oj ~ 
A described by the second term in the r.h.s. of Eq. (p^) . In particular, in the quasi- ID 
case the function /2 — /a determining the spatial dependence of this term has the form 



/2 — /s 



3g^ 
-^ - 2t> + x2 - 1/30. 



-'ir)--^(r 



(24) 



where B^{x) = x" 

Let us remind the reader that the above derivation is valid for u <^ Ec. In order 
to obtain the results in the range u > Ec one can calculate the sigma-model correlation 
functions entering Eqs. (0) by means of the perturbation theory 2°. We find then 



vmMr,)Mr2W)i 



V'{rkir^)M^,)Mr2)i^!ir2))i 



1 
+ 2 



Re{/i;d(r)n^(ri,r2) 



, (25) 



A;d(r) + Ren^(ri,r2), 



where n(^(ri,r2) is the finite-frequency diffusion propagator 

and the summation in Eq. (|26|) now includes q = 0. As was mentioned, the perturbation 
theory should give correctly the non-oscillatory (in a;) part of the correlation functions 
at a; ^ A. Indeed, it can be checked that Eqs. (p5D match the supersymmetric a- 
model results in this regime. Furthermore, in the 1/ g order [which means keeping only 



linear in n^^ terms in (p5D and neglecting —iuj in denominator of Eq.(p6D] Eqs.(p5D, 
( ^61) reproduce the exact results (piSf ), (|19|) even at small frequencies a; ~ A. We stress 
however that the perturbative calculation is not justified in this region and only the 
supersymmetry method provides a rigorous derivation of these results. 

As was mentioned, the case of broken time reversal symmetry (unitary ensem- 
ble) was considered above. Generalization to a system with unbroken time reversal 
symmetry is straightforward ^^; in l/g-oidei Eqs.(|13|), (H), and (p!9D are modified as 
follows: 

V^'(|^fe(ri)^fe(r2)|')ij = [1 + 2ka{r)] [1 + 2n^(ri, v^)] , (27) 

VWMvi)^i{v2)\^)e,u. - 1 = 2A;,(r)n^(ri,r2). (28) 

l^2(^*(ri)V'Kri)V'fc(r2)V'r (r2))i?,. = ^^(r) + [1 + k,{r)] Uoir^, ra), k^ I. (29) 

Using the supersymmetry method, one can calculate also higher order correla- 
tion functions of eigenfunction amplitudes. In particular, the correlation function 
(|V^^(ri)||'?/'|(r2)|)_B determines fiuctuations of the inverse participation ratio (IPR) 
I2 = /'^i'|'^^(r)|- Details of the corresponding calculation can be found in Ref.9; the 
result for the relative variance of IPR, 5{l2) = yaii^h) / {h)'^ being 

5{h) = -|/3 = 1^, (30) 

(3^ I3^g^ 

with a numerical coefficient ad depending on the sample dimensionality d and equal to 
ai = 1/90, 02 ~ 0.0266 and 03 ~ 0.0527 for quasi-lD, 2D, and 3D geometry respec- 



tively. The fiuctuations (30) have the same relative magnitude as the famous universal 
conductance fiuctuations. Note also that extrapolating Eq.(pO|) to the Anderson tran- 
sition point, where 5^ ~ 1, we find 5(/2) ~ 1, so that the magnitude of IPR fiuctuations 
is of the order of its mean value (which is, in turn, much larger than in the metallic 
regime; see the next section). 

Similarly, correlations of eigenfunction amplitudes determine fiuctuations of matrix 
elements of an operator of some (say. Coulomb) interaction computed on eigenfunctions 
ipk of the one-particle Hamiltonian in a random potential. Such a problem naturally 
arises, when one wishes to study the effect of interaction onto statistical properties of 
excitations in a mesoscopic sample (see below). 

Finally, the above consideration can be generalized to a ballistic chaotic system, by 
applying a recently developed ballistic generalization of the a-model ^^' ^^. The results 
are then expressed in terms of the (averaged over the direction of velocity) kernel 
5'(ri, ni; r2, n2) of the Liouville operator K = t>^nV governing the classical dynamics 
in the system, 

n(ri, r2) = J dnidn2 fi'(ri, rii; r2, 02); 

Kg{r:i_, rii; r2, n2) = (ttz/)"^ [(5(ri - r2)(5(ni - n2) - V"^] . (31) 
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Here n is a unit vector determining the direction of momentum, and normalization 
/(in = 1 is used. Equivalently, the function n(ri,r2) can be defined as 

n(ri,r2) = y rfty (ini^(ri,ni,t;r2) , (32) 

where g is determined by the evolution equation 

[^ + ^FniViU(ri,ni,t;r2)=0, t>0 (33) 

with the boundary condition 

~g\t=o={'^i^r'[S{r^-T2)-V-']. (34) 

Eq.(|3T|) is a natural "ballistic" counterpart of Eq.(|10|). In particular, generalization of 
Eqs.(|l3]), (|27| ) for the correlations of an eigenfunction amplitudes in two different points 
onto the ballistic case reads ^^ 

a(ri,r2,^) = l + -n(ri,r2). (35) 

Note that the ballistic cr-model approach is of semiclassical nature and thus valid on 
distances much larger than the wave length Xp. For this reason, Eqs.(^) represent 
the smoothed correlation function, which is not valid for |ri — r2| < Xp- Indeed, the 
contribution to n(ri,r2) from the straight line motion from r2 to ri can be easily 
evaluated, yielding e.g. in the 2D case n'^°)(ri,r2) = l/{'KpF\ri — r2|). This is nothing 
else but the smoothed version of the function kd{\r-i_ — r2\) = Jo(pF|ri — i"2|) giving 
the leading contribution to the short-scale correlations. A formula for the variance of 
matrix elements closely related to Eq.(^) was obtained in the semiclassical approach 
in Ref.25. In a very recent paper ^^ a similar generalization of the Berry formula for 
(V'fc(ri)'?/'fc(r2)) was proposed. 

Eq. (|35|) shows that correlations in eigenfunction amplitudes in remote points are 
determined by the classical dynamics in the system. It is closely related to the phe- 
nomenon of scarring of eigenfunctions by the classical orbits ^^' ^*. Indeed, if ri and r2 
belong to a short periodic orbit, the function n(ri,r2) is positive, so that the ampli- 
tudes |V'fc(ri)P and |'?/'i!c(r2)P are positively correlated. This is a reflection of the "scars" 
associated with this periodic orbits and a quantitative characterization of their strength 
in the coordinate space. Note that this effect gets smaller with increasing energy E of 
eigenfunctions. Indeed, for a strongly chaotic system and for |ri — r2| ~ L (L being 
the system size), we have in the 2D case n(ri,r2) ~ Xp/L, so that the magnitude 
of correlations decreases as E~^/'^. The function n(ri,r2) was explicitly calculated in 
Ref.l9 for a circular billiard with diffusive surface scattering. 



STRONG CORRELATIONS OF EIGENFUNCTIONS AND LEVEL RE- 
PULSION AT THE ANDERSON LOCALIZATION TRANSITION 

In the preceding section, we considered the metallic regime, where a typical wave- 
function V'i(r) is extended and covers uniformly all the sample volume. When system 
approaches the point of Anderson transition Ec-, these extended eigenfunctions become 
less and less homogeneous in space showing regions with larger and smaller amplitudes 
and eventually forming a multifractal structure in the vicinity of E^. To characterize the 



degree of non-homogeneity quantitatively, it is convenient to use the inverse participa- 
tion ratio (IPR) {I2{E)) = J dra{r,r,E) For extended states this quantity is inversely 
proportional to the system volume: {12(E)) = A{E)L~'^, with L and d standing for 
the system size and spatial dimension, respectively. The coefficient A in this relation 
measures a fraction of the system volume where eigenfunction is appreciably non-zero. 
In the regime of a good conductor A ^ 1 + 2/(3 ~ 1, whereas close to the mobility edge 
E = Ec it becomes large and diverges like A{E) (x \E — Ec\~'^^, with a critical index 
/X2 > ^^. This means that eigenfunctions become more and more sparse, when the 
system approaches the critical point of the Anderson transition. Just at the mobility 
edge eigenfunctions occupy a vanishing fraction of the system volume and IPR scales 
like {12(E)) oc L^"^^"^, with r/ > 0. Such a behavior reflects multifractal ^°' ^^' ^^ struc- 
ture of critical eigenstates. At last, in the insulating phase any eigenstate is localized 
in a domain of finite extension ^ and IPR remains finite in the limit of infinite system 
size L -^ 00. 

This transparent picture serves as a basis for qualitative understanding of spectral 
properties of disordered conductors. Indeed, as long as eigenstates are well extended and 
cover the whole sample, they overlap substantially and corresponding energy levels repel 
each other in the same way as in RMT. As a result, the Wigner-Dyson (WD) statistics 
describes well energy levels in a good metal ^^' ^°' ^^. In contrast, in the insulating 
phase different eigenfunctions corresponding to levels close in energy are localized far 
apart from one another and their overlap is negligible. This is the reason for absence 
of correlations of energy levels in this regime - the so-called Poisson statistics. 

However, a naive extrapolation of this argument to the vicinity of the transition 
point would lead to a wrong conclusion. Indeed, one might expect that sparse (mul- 
tifractal in the critical point) eigenstates fail to overlap, that would result in essential 
weakening of level correlations close to the mobility edge and vanishing level repulsion 
at E = Ec- However, numerical simulations show ^^' ^'^' ^^' ^^' ^^ that even at the mobil- 
ity edge levels repel each other strongly, though the whole statistics is different from 
the WD one. The purpose of this section is to explain how this apparent contradiction 
is resolved. We will show that critical eigenstates for nearby levels are so strongly 
correlated that they overlap well in spite of their sparse structure. 

To calculate the overlap function a{r,r,E,uj) in the critical regime ^^, we will 
exploit an exactly solvable model of the Anderson transition - so-called sparse random 
matrix (SRM) model ^^. This model is essentially equivalent to a tight-binding model, 
which is locally of tree-like stricture (i.e has no small-scale loops), with matrix size 
A^ playing a role of the system volume (number of sites of the tight-binding model). 
The SRM model can be used to construct an effective mean-field theory of Anderson 
localization ^^ corresponding to the limit d = 00. The inverse participation ratio 
{12(E)) = Na{r,r,E) is proportional to 1/A^ in the delocalized phase, as expected: 
{12(E)) = A{E)/N. The coefficient A{E) diverges close to the transition point, \E — 
Ec\ <^ Ec, like A(E) oc exp f constl-E — E^^^^'^j ^^ which differs from the power-law 
behavior expected for conventional rf— dimensional systems. The origin of such a critical 
dependence was explained in ^^' ^° and stems from the fact that A{E) is determined 
essentially by the "correlation volume" V{^) (i.e. number of sites at a distance smaller 
than correlation length ^), which is exponentially large, V{C,) oc exp(const^), for tree- 
like structures, whereas V{^) oc ^'^ for a d— dimensional lattice. Having this difference 
in mind, one can translate the results obtained in the framework of rf = 00 models to 
their finite-dimensional counterparts '^^. 

The calculation of the inverse participation ration Na{r,r,E) can be extended 
onto the overlap correlation function a{r,r,E,uj). Evaluating the r.h.s. of Eq.(^, we 
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find the following relation 



17. 



a(r, r, E, uj) = j—^a{r, r, E) . (36) 



This relation is valid everywhere in the phase of extended eigenstates, up to the mobility 
edge E = Ec, provided the number of sites (the system volume) exceeds the correlation 
volume. In particular, it is valid in the critical region \E — E^ ^ E^, where a typical 
eigenfunction is very sparse and a(r, r, E) grows like exp (const|ii^ — E'd^^^^j- 

Eq.(|36|) implies the following structure of eigenfunctions within an energy inter- 
val SE = UJ < A~^{E). Each eigenstate can be represented as a product \E'j(r) = 
^j(r)$£;(r). Here the function $£;(r) is an eigenfunction envelope of "bumps and dips" 
It is the same for all eigenstates around energy E, reflects underlying gross (multifrac- 
tal) spatial structure and governs the divergence of the inverse participation ratio (i.e. 
of the factor A{E)) at the critical point. In contrast, ipi{r) is Gaussian white-noise 
component fluctuating in space on the scale of lattice constant. It fills in the envelope 
function $£;(r) in an individual way for each eigenfunction, but is not critical, i.e. is 
not sensitive to the vicinity of the Anderson transition. These Gaussian fluctuations 
are responsible for the factor (3/ {(3 + 2) (which is the same as in the corresponding 
Gaussian Ensemble) in Eg . (|36|) . 

As was already mentioned, this picture is valid in the energy window 5E ~ A^^{E) 
around the energy E\ the number of levels in this window being large as 6E/A ~ 
NA~'^{E) ^ 1 in the thermodynamic limit A^ — > oo. These states form a kind of 
Gaussian Ensemble on a spatially non-uniform (multifractal for E -^ Ec) background 
<l>s(r). Since the eigenfunction correlations are described by the formula (^), which 
has exactly the same form as in the Gaussian Ensemble, it is not surprising that the 
level statistics has the WD form everywhere in the extended phase ^^. 

We believe on physical grounds that the same picture should hold for a conventional 
(i-dimensional conductor. First of all, the general mechanism of the transition is the 
same in d < cxd and d = oo models. Furthermore, the sparsity (multifractality) of 
eigenstates near the transition point takes its extreme form for d = oo models ^°, so 
that since the strong correlations (^) take place at c? = oo it would be very surprising 
if they do not hold at finite d as well. Finally, Eq.(|36|) is supported by the results of 
the calculations in the weak localization regime. Keeping only the leading (in d > 2) 
terms proportional to the powers of n(r, r) and considering the unitary ensemble for 
definiteness, we have up to the two-loop order (see Eq.(|^) and Ref.9) 



a{T,r,E,uj) = V~^ 



l + n(r,r) + ^n2(r,r) 



^a(r,r,E), (37) 



in full agreement with Eq. (p6|) . 

Replacing A{E) by the (i-dimensional correlation volume ~ C,'^, we conclude that 
for E close to Ec Eq.(PBD should be valid for u < A^, where A^ oc 1/^'^ is the level 
spacing in the correlation volume. For larger u, a{r,E,uj) is expected to decrease as 
ijj"'^!'^ according to the scaling arguments ^^' ^^' ■*^, so that we find cr(r, E^ u;)/a(r, E) ~ 
(a;/A^)~^/'^, up to a numerical coefficient of order of unity. Again, for any value of the 
energy E in the delocalized phase, taking the system size L large enough, L S> ^, we 
have a large number of levels 5E/ A ~ A^/A oc {L/^Y i^^ the energy window 6E where 
Eq . (|36|) holds, so that the level correlation will be of the WD form. 

Finally, let us consider what happens when we go from the critical regime (^ large, 
but L ::^ ^) to the critical point (^ :§> L). For this purpose, let us keep the system size 
L fixed and change the energy toward Ec^ so that ^ increases. When ^ gets comparable 
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to the system size, ^ ~ L, we have A^ ~ A. This is the border of apphcabihty of the 
above consideration. Correspondingly, we find 

or(r, r, E, uj)/a{r, r, E) ~ 1 , u < A (38) 

and a{r,r, E,uj)/a{r,r, E) ~ (cu/A)^''/'^ for tu > A. When E approaches further Ec, 
the correlation length C, ^ L gets irrelevant, so that these results will hold in the critical 
point (^ = oo). Of course, Eq.(PBD is not sufficient to ensure the WD statistics in the 
critical point, since there is only of order of one level within its validity range 6E ~ A. 
Indeed, the numerical simulations show that the level statistics on the mobility edge is 
different from the WD one ^^' ^'^' ^^' ^^' ^''. 

However, Eq.p8|) allows us to make an important conclusion concerning the be- 
havior of R2{uj) at small u < A, or, which is essentially the same, the behavior of 
the nearest neighbor spacing distribution P{s), s = cu/A, at s < 1. For this pur- 
pose, it is enough to consider only two neighboring levels. Let their energy differ- 
ence be ujq ~ A. Let us now perturb the system by a random potential V^(r) with 
(V(r)) = 0, {V{r)V{r')) = r6{r — r'). For the two-level system it reduces to a 2 x 2 
matrix {Vij}, i,j = 1,2, with elements Vij = /(i'^r V(r)\l'*(r)\E'j(r). The crucial point 
is that the variances of the diagonal and off-diagonal matrix elements are according to 
Eq. (|38|) equal to each other up to a factor of order of unity: 

{V,\)/{\V,\\) = air,r,E,uj)/air,r,E) ~ 1 (39) 

The distance between the perturbed levels is given by u; = [{Vn — V22 + (-i-'o)^ + |^i2p]"'^''^- 
Choosing the amplitude of the potential in such a way that the typical energy shift 
Vii ~ A and using Eq. (pQ]) , we find (|Vi2p) ~ A. As a result, the probability density 
for the level separation uj is for cj <^ A of the form dP ~ {uj/A)^duj/A, with some 
prefactor of order of unity. We thus conclude that in the critical point P{s) ~ c/ss^ 
for s <^ 1 with a coefficient cp of order of unity, in agreement with the numerical 
findings ^■*' ^^' ^''. 



FLUCTUATIONS IN THE ADDITION SPECTRA OF QUANTUM DOTS 

In this section we discuss effects of the eigenfunction statistics on the spectral 
properties of quantum dots. The electron levels of a quantum dot can be resolved if the 
temperature T is less than the mean single-particle level spacing in a dot and can be 
studied in transport experiments (see the recent review ^^ and references therein). In 
small dots containing only few electrons these levels show a regular structure familiar 
from the atomic physics ^^. What we have in mind here are however larger dots (con- 
taining in typical experiments from few hundred to few thousand conducting electrons). 
These dots are either disordered (diffusive) or, although being ballistic, are expected to 
have chaotic dynamics due to their irregular shape. As a result, it is natural to use a 
statistical description of the properties of energy levels and eigenf unctions in such dots. 

There are two types of the quantum dot spectra studied experimentally via mea- 
suring their I-V characteristics: (i) excitation spectrum, when excited levels are probed 
in a dot with given number of electrons by increasing the source-drain voltage, and (ii) 
addition spectrum, when electrons are added one by one by changing the gate volt- 
age. In the latter case, which will be the subject of our consideration here, one finds 
narrow conductance peaks separated the regions of (almost) zero current (Coulomb 
blockade). Statistics of the spacings between these peaks was studied in a number of 
recent experiments ^^' '^^' '^^; we will return to the experimental results below. 
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The simplest theoretical model which may be used to study distribution of the 
spacings is as follows. One considers a dot as a fixed size diffusive mesoscopic sample 
and assumes that changing a gate voltage by an amount 6Vg simply reduces to a uniform 
change of the potential inside the dot by a constant 'jSVg, with certain numerical 
coefficient 7 ("lever arm"). Such a model was used for numerical simulations of the 
addition spectra in Refs.'^'^' ^'^ . Below we consider the statistics of peak spacings within 
this model ^^, and later return to the approximations involved. We will neglect the 
spin degree of freedom of electrons of first; inclusion of the spin will be also discussed 
in the end of the section. 

The distance between the two consecutive conductance peaks is given by 

Sn = {En+2 — Em+i) — {Em+1 — En) 

= ^^%Xl-^^%^\ (40) 

where E^ is the ground state of a sample with A^ electrons. In the second line of 
Eq. ( ppp we rewrote S^ in terms of the Hartree-Fock single electron energy levels, with 
111 denoting the energy of the state jj^j in the dot containing i electrons. It is convenient 
to decompose Sn in the following way 

= E1 + E2 (41) 

The quantity E2 is the distance between the two levels of the same one-particle (Hartree- 
Fock) Hamiltonian Hn (describing a dot with A^ electrons) and is expected to obey 
RMT; in particular {E2) = A and r.m.s.(ii^2) = ^A with a numerical coefficient a of 
order of unity [a = 0.52 (0.42) for the orthogonal (resp. unitary) ensemble]. On the 
other hand, Ei is a shift of the level #(A^ + 2) due to the change of the Hamiltonian 
Hn — > Hn+i accompanying addition of the electron ^{N + 1) to the system. It can be 
in turn decomposed into the following three contributions ^^ 

El = e'/C + jdv {\i^l^,{v)\ + l^^+^Wl) 5U{v) 

+ Jdrdr'\^l^,{r)\\rN^,{r')\U.i\r-r'\) 
= E^ + e['^ + Ef ) (42) 

Here C is the dot capacitance, 6U{r) is the change of the self-consistent potential due 
to addition of one electron (i.e. difference in the self-consistent potential in the dots 
with N and N + 1 electrons), and [/^(r) is the screened Coulomb interaction (with the 
subscript k denoting the inverse screening length). In particular, in the experimentally 
most relevant 2D case (which we will consider below) and assuming a circular form of 
the dot with radius R, we have 

SU{r) = -^{R'-rr'^\ (43) 

while Uk is given in the Fourier space by f/K(q) = 27Te'^/e{q + k) with k = 27re^i^/e and e 
being the dielectric constant. The first term in Eq.(^2[) (the charging energy) determines 
the average value {Ei) and thus the average peak spacing (Sn) (since e^/C ^ A for a 
large dot with A^ ^ 1). This is the only contribution to Ei which is kept by so-called 
constant interaction model, which in addition neglects fluctuations of the capacitance 
C. Consequently, fluctuations of Sn in the constant interaction model are determined 
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solely by fluctuations of the single-particle level spacing E2 and thus should be described 
by RMT: r.m.s.(^7v) = aA. 

The term Ei in Eq. (|4l|) is however an additional source of fluctuations and is 
thus responsible for the enhancement of fluctuations in comparison with RMT. In 
principle, all three terms Ei , Ei , and Ei in Eq.(|^) contribute to this enhancement. 
Fluctuations of the flrst one, El = e^/C are due to the fact that the capacitance is 
slightly different from its purely geometric value because of a flnite value of the screening 
length. The corresponding correction to C can be expressed in terms of the polarization 
operator P(r, r') ^^. The latter is a fluctuating quantity (because of fluctuations of the 
eigenfunctions in the Fermi sea) and contains a random part Pr{r,r') leading to the 
following expression for the random part of the charging energy: 



(eVC), = 2 f drdr'6U{r)Pr{r,r')6U{r'). 



Evaluating the fluctuations of the polarization operator ^^, we flnd f\ 



(44) 



vaT{E[ 



(0)^ 




n2 



rfrirfr25t/(ri)n(ri,r2)<5f/(r2) 



(45) 



Now we consider fluctuations of the last term. El , in Eq. (^2]) . Using Eqs.(|13[), 
for the correlations of eigenfunction amplitudes in two remote points, the variance 
of El is found to be 

var(sP) = -i;;j|dridr;dr2rfr'2f/«(|ri-r;|)?7«(|r2-r^|)n(ri,r2)n(r;,r^) 



oc 



4A^ 



P'\9, 



J dr,dr.U^r,,r.) 



(46) 



Finally, fluctuations of the term Ei can be also evaluated with help of Eqs.(|T3D, (pTf), 
yielding 



var(£'} 



(l)^ 



OC 



1A2 



rfridr25f/(ri)n(ri,r2)5f/(r2) 



(47) 



It is seen that ioi g ^ 1 all the contributions Eqs.(^5D-(^7D are parametrically 
small compared to the RMT fluctuations (which are ~ A). Fluctuations of the term 
E[ related to the change SU{r) of the self-consistent potential represent parametrically 
leading contribution to the enhancement of the peak spacing fluctuations with respect 
to RMT. 

Let us now discuss approximations made in the course of the above derivation: 



''We use the obvious notations for the variance and the root mean square deviations of a quantity X: 
var(X) = (X^) - (X)2; r.m.s.(X) = [var(X)]i/2. 



12 



i) It was assumed that changing of the gate voltage results in a spatially uniform 
change of the potential in the sample, that led us to the expression Eg . (^31) 
for the change of the self-consistent potential 6U{r) accompanying the addition 
of one electron to the dot. This result would correspond to a gate located far 
enough from the sample. In a more realistic situation, when the gate is relatively 
narrow and located close to the sample, the potential change SU{r) (as well as the 
additional electron density) will be mainly located on the side of the dot facing 
the gate. Furthermore, the total area of the dot is not fixed, so that the dot 
gets larger (and slightly deformed) with adding each electron to it. These effects 
lead to some increase of the fluctuations of El in comparison with the model 
considered above. If the size of the gate and its distance to the dot are of the 
same order of magnitude as the dot size, then 

(1)^ 



T.m.s.{El') oc A/ 

as in Eg . (^T]) , with a geometry dependent-numerical prefactor. The upper bound 
for the magnitude of fluctuations of E^ is given by the opposite limiting case, 
when additional electron density occupies an area ^ Xp x Xp near the gate, in 
which case one finds 

T.m.s.{E[^^) ~ A/y^. 

ii) The dot was supposed to be diffusive in the calculation. For a ballistic dot one 
should replace n(r, r') by IlB{r,r'), as was explained above. This would mean 
that the parameter g is replaced by ~ A^^/^ ~ L/Xp, where A^ is the number 
of electrons in the dot and L the characteristic linear dimension. The numerical 
coefficient would depend, however, on "how strongly chaotic" is the dot. Role of 
the eigenfunctions fluctuations and correlations ( "scars" ) in enhancement of the 
peak spacing fluctuations was studied in ^^ via numerical simulations of a dot 
with A^ ^ 100 electrons. 

iii) It was assumed that the dot energy and the measured gate voltage are related 
through a constant (or smoothly varying) coefficient 7. This "lever arm" 7 de- 
pends, however on the dot-gate capacitance, which is also a fluctuating guantity. 
If the gate size and the distance to the gate is of the same order as the size of 
the dot, these fluctuations should be of the same order as fluctuations of the dot 
self-capacitance given by Eg. (^5[) , and thus lead to additional fluctuations which 
are parametrically small compared to A. In a more general situation (thin gate 
located close to the sample) an additional analysis along the lines of Ref.l6 is 
necessary. 

iv) The calculation was done within the random phase approximation, which assumes 
that Ts = e^/evp ^ 1, with vp being the Fermi velocity. In realistic dots however 
Ts ~ 1. Since this value is still considerably lower than the Wigner crystallization 
threshold, the calculations should be still valid, up to a numerical factor a{rs) 
[depending on r^ only and such that a{rs -C 1) = 1]. 

v) We considered the model of spinless electrons up to now. Let us briefly discuss 
the role of the spin degree of freedom. Within the constant interaction model, it 
would lead to a bimodal distribution ^"^ of peak spacings 



nsN) = I 



S{Sn - e'/C) 




(48) 
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where Pwd{s) is the Wigner-Dyson distribution and A denotes the level spacing 
in the absence of spin degeneracy. The value of the coefficient a in the relation 
r.m.s.(5'Ar) = aA is then increased (compared to the spinless case) and is equal 
to 1.24 (1.16) for the orthogonal (resp. unitary) ensemble. Taking into account 
fluctuations of eigenfunctions (and thus of Ei) however modifies the form of the 
distribution ^^. The value of the term i^j representing the interaction between 
two electrons is larger in the case when ipN+2 and ipN+i correspond to two spin- 
degenerate states (i.e. have the same spatial dependence of the wave function), 
since 



drdr'\^Ur)\\^f{T')\U^{\T-r'\)) - {j drrfr'|^f (r)||^|(r')|f/«(|r - r'|)) 

J drdr'kd{\r - r'|)f/«(|r - r'|) ~ A (49) 



for Ts ~ 1 (the coefficient depends on r^, see ^^). Therefore, filling a state ipi^ 
pushes up the level ipii (with respect to other eigenstates) by an amount of order 
of A. This removes a bimodal structure of the distribution of peak spacings and 
slightly modifies the value of the coefficient a. 

Basing on the above analysis, we can make the following general statement. Imag- 
ing that we fix r^ ~ 1 (i.e. fix the electron density and thus the Fermi wave length) and 
the system geometry, and then start to increase the linear dimension L of the system. 
Then, while the average value of the peak spacing S^ scales as (S^) ~ e^/C oc 1/L, 
its fiuctuations will scale differently: r.m.s.(S'7v) ~ A oc 1/L^. This result is not at all 
trivial, since in an analogous problem for classical particles '^^' ^° the fiuctuations are 
proportional to the mean value {3^). The physical reason for smaller fiuctuations in 
the quantum case is in the delocalized nature of the electronic wave functions, which 
are spread roughly uniformly over the system. 

The above prediction was confirmed by a recent experiment ^^, where a thorough 
study of the peak spacing spacing statistics was carried out. It was found that the low- 
temperature value of r.m.s.(S'Ar), as well the typical temperature scale for its change are 
approximately given by the mean level spacing A (while in units of Ec the magnitude 
of fiuctuations was very small, typically 2-4%). We note also that in recent numerical 
simulations ^^ fiuctuations of the addition energies were found to be approximately 
0.7A, in agreement with our results. 
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